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Abstract 

We present a strategy designed to separate several possible origins of the well-known 
enhancement of the AI = 1/2 amplitude in non-leptonic kaon decays. In particular, 
we seek to unambiguously quantify the role of the charm quark mass in the observed 
enhancement. This is achieved by considering QCD with an unphysically light charm 
quark, so that the theory possesses an approximate SU(4)l x SU(4)r chiral symmetry. 
The strategy proceeds by computing the relevant operator matrix elements and monitor- 
ing their values as the charm quark mass departs from the SU(4)-symmetric situation. 
We study the influence of the charm quark mass in Chiral Perturbation Theory. First 
results from lattice simulations in the SU(4)-symmetric limit are also discussed. 
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1 Introduction 

A quantitative understanding of non-leptonic kaon decays, such as K — > tttt, has been 
elusive for many years, and thus the explanation of the famous AI =1/2 rule or the 
value of e'/e has remained a longstanding problem. Early analyses have shown that, 
in a Standard Model-based explanation, the bulk of the enhancement of the AI = 1/2 
amplitude must be due to long-distance contributions generated by the strong inter- 
actions [1]. A reliable determination of these effects must inevitably be performed at 
the non-perturbative level [2,3]. Currently, lattice simulations of QCD are the only 
known methodology which can achieve this goal with controlled systematic errors. Lat- 
tice studies of K — > tttt decays have, however, been hampered by a number of technical 
difficulties, for instance by the so-called Maiani- Testa No-Go Theorem [4], which states 
that transition amplitudes of two-body decays cannot simply be obtained from the 
asymptotic behaviour of Euclidean correlation functions computed in very large vol- 
umes. In addition, lattice calculations employing commonly used Wilson fermions need 
non-perturbative subtractions of power divergences for defining properly renormalised 
operators which enter the effective weak Hamiltonian [5]. 1 Thus, efforts to investigate 
the AI =1/2 rule numerically on the lattice had practically come to a halt during most 
part of the 1990s. 

Recently, it was realised, though, that the Maiani- Testa Theorem does not apply 
to volumes with linear extensions of a few fm, and a relation between the physical two- 
body decay rate and the square of the transition matrix element in finite volume could 
be derived [8] (for subsequent work, see [9]). Furthermore, the advent of fermionic 
discretisations which preserve chiral symmetry at non-zero lattice spacing (Ginsparg- 
Wilson fermions) [10-18], has alleviated the mixing problem to the extent that the 
renormalisation patterns of four-fermion operators mediating AS = 1 transitions are 
like in the continuum [19]. These developments have removed the main obstacles for a 
lattice treatment of non-leptonic kaon decays, based on first principles. An approximate 
realisation of Ginsparg- Wilson fermions via the so-called Domain Wall formulation has 
already been applied to non-leptonic kaon decays with some success [20,21]. 

The decay of a neutral kaon into a pair of pions in a state with isospin / is described 
by the transition amplitude 

T(K° - H I=a ) = iA a e iS « , a = 0, 2, (1.1) 

where 8 a is the scattering phase shift. In this paper we focus on the AI =1/2 rule, i.e. 
the empirical observation that the amplitude \Aq\ is significantly larger than | ^4.2 1 ? 

|4)|/|42| = 22.1. (1.2) 

There are several possible sources for an enhancement due to strong interaction effects. 
These involve physics at the scale of the charm quark, i.e. at around 1 GeV, physics at an 

1 For recent progress using twisted mass QCD, see [6] and [7]. 
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intrinsic QCD scale Eqqb ~ 250 MeV, and effects due to pionic final state interactions 
at around f 00 MeV. It remains unclear, though, whether the enhancement is the result 
of an accumulation of several effects, each giving a moderate contribution, or whether 
it is mainly due to a single cause or mechanism. 

The possible role of the charm quark has been pointed out long ago in ref. [22]: 
for energy scales E <C m c the GIM mechanism no longer operates, which gives rise 
to so-called "penguin operators" that mediate AI = 1/2 transitions. The commonly 
accepted scenario that the AI = 1/2 rule arises predominantly from long-distance QCD 
contributions has been supported by several analytical [23-27] and numerical [20, 21] 
studies. 

In this paper we describe a general strategy which allows to disentangle contri- 
butions from the various sources and quantify them using numerical simulations. The 
main idea is to keep an active charm and determine the leading low-energy constants 
(LECs) associated with the CP conserving AS = 1 weak Hamiltonian of the chiral low- 
energy effective theory as a function of the charm quark mass. These parameters can be 
determined by computing suitable correlation functions at small masses and momenta 
in a lattice simulation and comparing them with the predictions of Chiral Perturbation 
Theory (ChPT). If the charm is degenerate with the light quarks, m u = m<j = m s = m c , 
the theory has an exact SU(4)l X SU(4)r symmetry in the chiral limit, which is bro- 
ken explicitly by the weak interactions. In this situation the calculation of the LECs 
corresponding to AI =1/2 and 3/2 transitions will expose any QCD contribution to 
the enhancement at scales around -Eqcd- The effect of a heavier charm quark can 
then be isolated by monitoring the amplitudes as m c departs from the degenerate limit: 
m c > m u = nrid = m s . By keeping an active charm quark and employing a lattice for- 
mulation which preserves chiral symmetry, power-divergent subtractions of the relevant 
operators can be avoided, and small quark masses can be simulated without numerical 
instabilities. 

Since simulations with dynamical Ginsparg- Wilson fermions are still prohibitively 
expensive, it is reasonable to perform initial tests of our proposed strategy in the 
quenched approximation. It is well known, though, that the quenched theory is af- 
flicted with several problems, whose implications are discussed in detail in the relevant 
sections below. In order to illustrate our strategy, we present first numerical results 
for the SU(4)-symmetric case obtained in the quenched approximation for quark masses 
near m s /2. Numerical results for smaller masses and the investigation of the dependence 
on the charm quark mass are left to future publications. 

The outline of the remainder of this paper is as follows: in Sect. El we discuss the 
effective weak interactions with an active charm quark. The lattice-regularised theory, 
formulated using Ginsparg- Wilson quarks, is described in Sect. EH The renormalisation 
group invariant formulation of the effective weak Hamiltonian is addressed in Sect. 0J 
In Sect. El we discuss the effective low-energy description of AS = 1 weak interactions 
in terms of ChPT in a finite volume. The decoupling of the charm quark for small 
quark masses, i.e. m c > m u ,m^,m s , is analysed in Sect. Housing ChPT. Our numerical 
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results are discussed in Sect. [JJ and in Sect. |H] we present our conclusions. In several 
appendices we provide further information on the SU(4) classification of operators, the 
transformation properties of four-quark operators in the lattice theory, the perturbative 
renormalisation of four-quark operators using overlap fermions, as well as details of the 
calculations performed in finite-volume ChPT. 

2 The effective weak interactions with an active charm quark 

In this section we collect some basic facts and definitions regarding the AS = 1 ef- 
fective weak interaction, focusing on the less familiar case of an active charm quark. 
Throughout this paper we work in Euclidean space-time. 

2.1 Operator product expansion and global symmetries 

The decay of a (neutral or charged) kaon into two pions is induced by charged-current 
weak interactions, mediated via the exchange of a VF-boson. It can be described in 
terms of an effective V — A current-current interaction, i.e. 

= \gl £ (V qs )*V qd J d 4 xd 4 y (s^P-q)(x) D^{x - y) (q lu P-d)(y), (2.1) 

q=u,c,t 

where g\ = 4v / 2GfM^ / , V qs , V q d denote elements of the CKM matrix, P_ = ^(1 — 75), 
and D^ v is the propagator of the VT-boson. Contributions from the top-quark are 
suppressed by three orders of magnitude relative to those from the up-quark and can 
thus be safely neglected. At this level of accuracy one has (V us )*V u d = — (V cs )*V C( i, so 
that 

S w = lgl(V us )*V ud J d 4 xd 4 y 

x{(s^P-u)(x)D^(x - y) (urf v P-d)(y) -(u^c)}. (2.2) 

The dominant contribution to S w comes from the region x ~ y, so that the integral can 
be evaluated using the operator product expansion: 

S w « / d 4 xH w (x), H w (x) = -&(v us yv ud Y,KQn{x), (2.3) 

where the coefficients k n depend on the VK-boson mass, My/, and the renormalisation 
scheme used to define Q n . In order to classify the operators that can occur in the sum, 
we now discuss the global symmetries that must be respected. 

To this end we consider QCD for two generations and write its action as 

S = S G + J d 4 x^(x) (D + MP + + M f P_) ip(x), (2.4) 

where Sq denotes the gauge action, D the massless Dirac operator, M the quark mass 
matrix, and the four-flavour quark field, with flavour components u, d, s, c. The action 
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is invariant under SU(4)l x SU(4)r chiral transformations if M is transformed according 
to the 4 ® 4* representation. For real diagonal M it is also invariant under parity and 
charge conjugation. 

Under SU(4)l x SU(4)r the effective action 5 W divides into two parts that belong 
to different representations. The decomposition reads 

q _ q+ I C- 

= \gl{V U s)*V ud J d A xd 4 y{[V suud (x,y) ±V sudu (x,y)} - (u -> c)}, (2.5) 

where 

Pafasix, V) = {^alfiP-^)^) D^(x - y) (ippj v P-^s)(y), (2.6) 

with generic flavour indices a, ... ,5. From the structure of V a p^s it is evident that 
both parts are singlets under SU(4)r, while it can be shown that and S~ trans- 
form according to the irreducible representations of SU(4)l of dimensions 84 and 20 
respectively [28]. Furthermore, the effective action, as well as the QCD action, are in- 
variant under the CPS transformation, i.e. combined operations of CP followed by an 
interchange of d and s quarks [3] . 

Having thus listed the global symmetries, we now proceed to find the operators 
of dimension d < 6 which occur in the operator product expansion, eq. (j2,3[) . noting 
that operators of higher dimension are suppressed by powers of M\y. We start the 
discussion with four-quark operators; given a suitable basis of operators, we seek to 
construct linear combinations that are SU(4)r singlets and which transform according 
to irreducible representations of SU(4)l- The result is 

Qf = {[OiUud ± PiUdu) - (u ->• c), (2.7) 
[OiW = &alnP-^)(^p^P-^s). (2.8) 

Note that the decomposition yields unique operators Qf , Qj - , which transform accord- 
ing to irreducible representations of dimensions 84 and 20, respectively. Further details 
are described in Appendix [X] Furthermore, the requirement of CPS symmetry is also 
fulfilled. 

Two-quark operators of dimension d < 6 which are singlets under SU(4)r can be 
constructed by including derivatives and two or more factors of the mass matrix. Drop- 
ping operators that are related via the equations of motion and operators transforming 
according to low-dimensional representations of SU(4)l one is left with [30] 

{MM\0 a P + [Mi>} & ) and (MM%^M^] a P^ 6 ). (2.9) 

The projected operators that transform according to representations of dimension 84 
and 20 and which satisfy CPS symmetry are then 

Qf = {[0 2 ]suud ± [02]sudu) ~ (u c), (2.10) 

[02^6 = (MM%{® a P + [Mi)]s) + ([?Aft] Q P_^)} 

+ (<*«-► /J, 7 «->$). (2.11) 
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For M = diag(m u , raj, m s , m c ) the expression for Q 2 is 

Qt = (m 2 u - m 2 c ){m d (sP+d) + m s (sP^d)}. (2.12) 

Zero-quark operators which transform under SU(4)l can also be constructed using at 
least four powers of the mass matrix. Any such terms conserve flavour automatically 
(the mass matrix is diagonal) and hence they do not contribute to the part of the 
effective weak Hamiltonian we are interested in. 

2.2 Renormalisation and mixing 

So far we have neglected the issue of renormalisation. We note that the effective action 
S w in eq. (|2.1j) is finite and does not need any subtractions. This is obvious from the 
fact that the charged currents are not anomalous in QCD and thus have a natural 
normalisation. After the operator product expansion, the operators Qf and Q 2 of 
dimension 6 that appear in eq. ()2.3|) must be renormalised, and there is no reason why 
operators with the same transformation properties should not mix. Assuming that the 
adopted regularisation preserves enough of the relevant symmetry structure to exclude 
any other mixings, the relations between renormalised and bare operators are of the 
general form 

n ± _ 7 ± n ±,bare . j,± n ±,bare 

Qf = Z±Qf> haie + Zt 2 Qt' haxe . (2.13) 

The lattice formulation with Ginsparg- Wilson fermions discussed in the next section is 
an example of a regularisation where these relations hold without modification. 

Since the operators Q 2 ' are , defined as in eq. (|2,12|) . are linear combinations of the 
non-singlet chiral densities, which are multiplicatively renormalisable, we are free to 
set Z 2l = 0. The renormalised operators Q 2 are then obtained as linear combinations 
of renormalised densities with coefficients that are polynomials in the renormalised 
quark masses. As is obvious from eq. (|2.12j) . the Glashow-Iliopoulos-Maiani (GIM) 
mechanism ensures that any contribution proportional to Q 2 vanishes when m c = m u . 
In particular, the mixing of Q 2 with Qf is absent in the chiral limit, and hence the 
factors can be fixed at vanishing quark masses. The constants Z^ 2 must then be 
determined by requiring that any residual divergences in the matrix elements of for 
m u 7^ m c are cancelled. 

We note that the mixing of Q 1 with Q 2 is usually ignored in matrix elements in 
which there is no momentum transfer between initial and final states, \i) and |/). The 
reason for this is that the operators (~sP±d) can be written as a total four-divergence, 
and hence the matrix element 

(f\(sP±d)\i) = 

(m s - ma)- 1 d^(f\(s^d)\i) ± (m s + m d )~ l d^f^^d^i) (2.14) 
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vanishes if \i) and |/) have the same four-momentum. In other words, only the terms 
proportional to Q^ 1 contribute in physical matrix elements. Evidently this argumenta- 
tion is correct only if m s ± md / and if kinematical singularities at zero momentum 
transfer can be excluded. 



3 The effective weak Hamiltonian in lattice QCD 

The construction of the effective weak Hamiltonian on the lattice proceeds by finding a 
linear combination of local composite fields with coefficients such that the correct opera- 
tor (including normalisation) is obtained in the continuum limit. If a lattice formulation 
is used which preserves chiral symmetry, many of the difficulties encountered with Wil- 
son fermions, such as the mixing with lower dimensional operators, can be avoided [19]. 
In particular, one may require that the operator basis transforms in a simple way under 
the chiral SU(4)l x SU(4)r symmetry group at non-zero lattice spacing. 

3.1 Ginsparg- Wilson fermions 

The formulation of lattice QCD with exact chiral symmetry proceeds by introducing a 
lattice Dirac operator D which satisfies the Ginsparg- Wilson relation 

75D + D 75 = aDj 5 D. (3.1) 

Explicitly we take D to be the Neuberger-Dirac operator [16], defined by 

D = -\l-A(AU)- 1 / 2 }, A = l + s-aD w , (3.2) 

where D w denotes the massless Wilson-Dirac operator, a is the lattice spacing, and s 
is a free parameter in the range \s\ < 1, which can be tuned to optimise the locality 
properties of D [18]. If a is set to 

S = TT7 (3 ' 3) 

it is straightforward to check that D satisfies the Ginsparg- Wilson relation. The in- 
finitesimal chiral transformations of the quark fields tp and Y> are given by [17] 

Si/j = A a 75 (l - aD)ifj, 5$ = ^ 75 A a , (3.4) 

where A a is a flavour matrix. Furthermore, it is useful to define the modified fermion 
field ip by 

$=(\-\aD)ij), (3.5) 
whose infinitesimal chiral transformation is given by 

5$ = A a 75^. (3.6) 

An important reason for introducing the modified field ip is that local composite oper- 
ators constructed from ip and ip have simple transformation properties under the chiral 
symmetry, similar to those in the continuum. 
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3.2 Operator basis 

In the construction of the weak Hamiltonian in the lattice theory we are concerned 
with finding a basis of local operators of dimension d < 6 with the same transforma- 
tion behaviour as Q\ and Q2 under the exact symmetries of the lattice theory. These 
include the discrete lattice symmetries, the gauge transformations, SU(4)l x SU(4)r 
chiral transformations, and the CPS symmetry. 

Concentrating first on four-quark operators, we note that in the continuum the 
problem is solved by constructing a basis of gauge-invariant operators that transform 
as scalar fields under the (restricted) Lorentz group SO(4). For Wilson-type fermions 
the classification proceeds along the same lines, except that the Lorentz group SO (4) 
is reduced to the hypercubic group SO(4, Z). More precisely, the task is to find all 
tensors that are invariant under the spin covering Spin(4, Z). A detailed analysis [31] 
then shows that no additional invariant tensors can occur under Spin(4, Z), and hence 
the basis of four-quark operators in the continuum theory is also a basis in the lattice 
theory. 

It is easy to convince oneself that four-quark operators on the lattice composed in 
terms of if] and the modified fields if) have exactly the same transformation behaviour 
under SU(4)l x SU(4)r as the corresponding operators in the continuum. Thus, the 
lattice counterpart of the operator Qf in eq. Q2.7JI can be chosen as 

Q ±,barc = ^^p_u)(u^ P _d) ± (s 7M P_d)(^P_£t)} - («-> c). (3.7) 

In order to find the set of two- and zero-quark operators in the lattice theory one has 
to classify the appropriate tensors according to their transformation properties under 
the spin covering of the hypercubic group. Following similar arguments as in the case 
of four-quark operators, one is left with only one operator, namely 

Q±' bare = (ml - m 2 c ){m d (sP + d) + m a (aP_d)}, (3.8) 

where m u , . . . , m c are the bare masses that appear in the lattice action. 

It has been pointed out that the infinitesimal chiral transformations in eq. Q3.4|) 
do not commute with CP [32]. This would imply that the operators in the basis have 
simple transformation properties either under the chiral symmetry or under CP, but not 
under both symmetries simultaneously. The discussion in [31], which is summarised in 
Appendix[Bl shows, however, that simple CP transformation properties are recovered if 
one considers insertions of Q 1 ' are and Q 2 ' 3X6 i n correlation functions of local operators 
at non-zero distances. In this situation, the operators transform under CP like in the 
continuum theory, up to an overall factor that depends on the bare quark masses. This 
is perfectly adequate for the study of issues like operator mixing. 

The upshot of this discussion is that the use of Ginsparg- Wilson fermions yields 
an operator basis in the lattice-regularised theory, whose mixing patterns are exactly 
like those found in the continuum. In particular, the GIM cancellation of contributions 
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proportional to two-quark operators is quadratic in the masses, and thus the coefficients 
that quantify the mixing of Q^' bare with (sP±d) cannot develop any power divergences, 
which are so hard to control for ordinary Wilson fermions. 

The derivation of the effective weak Hamiltonian and the discussion of operator 
mixing in the continuum (see Sect. I^J) assumes that a regularisation which preserves 
chiral symmetry exists. Lattice QCD with Ginsparg- Wilson fermions is, in fact, the 
only known such regularisation, so that the above discussion establishes the findings in 
the continuum theory in a rigorous manner. 

4 Renormalisation group invariant formulation 

We are now in a position to specify the effective weak Hamiltonian, which is given as 

n - = Sr^y^i E W Si + k 2 Q2} ■ (4-1) 

Here it is assumed that the operators Qf, are renormalised in a particular scheme 
and at a given value of the renormalisation scale. The dependence on the scheme and 
scale can be eliminated by passing to the so-called renormalisation group independent 
(RGI) normalisation, for which correlation functions of the operators stay unchanged 
along the renormalised trajectory. This requires knowledge of the anomalous dimen- 
sions, which, in the case of Qf, have been determined to two loops in perturbation theory 
for schemes like dimensional reduction (DRED) [33], 't Hooft-Veltman (HV) [34], naive 
dimensional regularisation (NDR) [34], as well as the regularisation-independent (RI) 
scheme [35,36]. To be more explicit, consider Qf in some renormalisation scheme (e.g. 
RI for definiteness) at scale n. The RGI counterpart of Qf is obtained via 

(Qf)R G i = c ± ^/A)Qf^), (4.2) 

where 

c+Ou/A) = (2b g 2 (riy° /{2bQ) exp j - ^ dg 

and g([i) denotes the QCD running coupling constant at scale /i. Here 7 ± (s0 and f3{g) 
are the anomalous dimension of Qf and the RG /3-function in the chosen scheme, with 
the respective one- loop coefficients 7^ and 60 - 2 Some perturbative results for the one- 
and two-loop coefficients are summarised in appendix O 

The renormalisation of four-quark operators for lattice actions which preserve chiral 
symmetry has been studied by a number of authors, both in perturbation theory [19,38- 
40], and also non-perturbatively [21,41]. For example, in ref. [19] the renormalisation 
factors Zf' , which relate matrix elements of Q^' bare in the RI scheme to those computed 

2 There is no standard convention for the overall normalisation of RGI operators in the literature. 
The one adopted here is similar to that used in [37] but differs, for instance, from the commonly used 
normalisation for the RGI kaon B-parameter. 



0{g) ' b g 



(4.3) 
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using overlap fermions, were determined in one-loop lattice perturbation theory. The 
RGI matrix elements can then be obtained via the anomalous dimensions, which are 
known to two loops [35,36]. 

The coefficients k^ 1 , k^ can be considered for the RGI operators or for those in a 
given continuum scheme. In either case the coefficients kf have a well-defined perturba- 
tive expansion. The coefficients k^ 1 , though, remain undetermined, because the effective 
Hamiltonian is derived from the fundamental theory by matching a set of on-shell am- 
plitudes, in which the total momentum that flows into the interaction vertex vanishes, 
while k,2 multiplies a term which can be written as a total four-divergence. It should 
be noted, however, that does not contribute to the physical K — > irir amplitude. 

If one adopts the RGI scheme for according to eq. (|4.2|) then the perturbative 
expression for kf 1 is obtained as 



kl{M w ) = (2b g\M w ))-^/^o) h + k <y)g^ Mw ) + 0(g 



a 



(4.4) 



If one follows the conventions of the MS-scheme to define the coupling g, then the 
expression for the coefficient k a] ^ reads 



^;(i) 



(4tt) 2 AN C (11N C _ 2jV f ) 2 

1229iV 2 + 1881 - N t ( 418iV c + 



693iV c + 126iV f 



-a 



450 \ 
~Nc) 



+ 80iY f 2 



(4.5) 
N c , and the 



where we have made explicit the dependence on the number of colours, 
number of active quark flavours, Nf. 

In the literature one often finds discussions of the short-distance contribution to 
the AI = 1/2 enhancement. Since Q± mediates AI =1/2 transitions exclusively, one 
considers the ratio 



R(jm) 



k^(M w )kf(p) 

kt(M W )k^(iJ,y 



m c < n < Mw- 



(4.6) 



The effects of physics above the charm quark mass can be estimated by computing 
R{m c ) via an integration of the perturbative /3-function for the coupling. One then 
obtains R(m c ) ~ 2, which is often regarded as a "first step" in the explanation of the 
AI = 1/2 rule. It must be kept in mind, however, that this analysis is not on a solid 
footing, since it relies on the applicability of two-loop perturbation theory down to scales 
where corrections are of order 100%. 



5 Weak interactions in finite-volume Chiral Perturbation Theory 

We now turn to the discussion of the A5 = 1 weak interactions in the framework of 
Chiral Perturbation Theory (ChPT). In particular, we will list expressions of correlation 
functions involving the counterparts of the operators Qf , in the effective low-energy 
description. 



9 



5.1 Weak Hamiltonian in the SU(4) chiral effective theory 

If one considers the unphysical situation of a light charm quark, QCD can be described at 
low energies by an effective Lagrangian which possesses an SU(4)l x SU(4)r symmetry. 
At leading order it is given by 

£ E = \F 2 Tr [{d^U)d^J ] } - |S Tr [uM^e id/N{ + Mt/V^l , (5.1) 

where U £ SU(4) denotes the field of Goldstone bosons, 9 is the vacuum angle, and 
M is the quark mass matrix. Although we are dealing with the SU(4)-symmetric case, 
we have explicitly indicated the iVf-dependence in the phase factor. At this order two 
effective coupling constants ("low-energy constants", LECs), F and S, appear, which 
denote the pion decay constant and the chiral condensate in the chiral limit, respectively. 
In order to incorporate the weak interactions, we need to find low-energy transcriptions 
of the operators which appear in the weak Hamiltonian, eq. (|4.1|) . in terms of the field U. 
To this end we introduce the left-handed current in the low-energy theory as 

J a ^\F\T a U(ud^) pa , (5.2) 

where the matrices T a denote the (hermitian) generators of SU(4)-flavour. The current 
J® is formally obtained by promoting the partial derivative in eq. (|5.1[) to a covariant 
one via 

d„U — D^U = (fy + iA°T a )U, (5.3) 
where is an external, left-handed flavour gauge field, such that 3 

By following a similar procedure in the fundamental theory it is easy to see that is 
the low-energy counterpart of 

J^ = ^P-(T a )a^. (5.5) 

Furthermore, by writing the expression for the four-quark operator [Oi] a /3 7 5 in eq. Q2.8|) 
in terms of products of currents J®, we can determine its representation in the low- 
energy theory (as far as symmetries are concerned) as 

[Oi W = i fA { Ud ^) ia ( ud ^) S p ■ (5 - 6) 

At leading order in the chiral expansion one can show that this is the only operator 
with the same symmetry properties as its counterpart in the full theory. In a similar 
manner one obtains a representation of the two-quark operator [O2] as 



[0 2 W = -|S(MMt) 7/3 ([7Mte i0 / JV f + MU^e~ i0 ^ 



6a 



+ (a^/3, 7 ^5). (5.7) 



Note that with this definition the current is formally imaginary. 
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The projection of 0\ and O2 onto operators 0^ and 2 , which transform under irre- 
ducible representations of dimensions 84 and 20, is performed according to the procedure 
outlined in Appendix The effective weak Hamiltonian at leading order in ChPT is 
then given by 

W^ hPT = J 2 (V us )*V ud ^ \ 9l C a p-jS{Ol} a/3-yS + 9*2 c a/3 7 <5 [Of ] a j3yS \ , (5-8) 
Z1V1 W a =± 

where g^, g 2 are low-energy constants and the coefficients c Q , / 3 7 ^ are Clebsch-Gordan- 
type numbers. For the physical flavour assignments chosen in eqs. ()2.7fl and (|2.1U|) . their 
values are given by 

Csuud — 1) C-sccd — 1) (5-9) 

with all other coefficients set to zero. 

At leading order in ChPT the ratio of amplitudes corresponding to AI = 1/2 and 
3/2-transitions is given by 

4> + (5,io) 

A non-perturbative determination of the ratio /gf thus yields direct information on 
the AI =1/2 enhancement. In the following we describe how the low-energy constants 
gf 1 and g 2 can be computed by matching suitable correlation functions evaluated in 
QCD to their analytic expressions obtained in ChPT. 

5.2 Chiral Perturbation Theory in finite volume 

Our task of determining the low-energy constants in eq. (|5.8|) through a comparison of 
correlators evaluated in lattice QCD with the expressions of ChPT can only succeed if 
the latter are valid in the parameter range accessible by lattice simulations. Numerical 
simulations of lattice QCD are necessarily performed in a finite volume, and therefore 
the approach to the chiral limit, where the predictions of ChPT are most accurate, 
will inevitably lead to strong finite-volume effects. However, within the framework of 
ChPT it is possible to account for finite-size effects in a systematic manner. Although 
our actual simulations, which we describe in Sect. do not yet reach the limit of very 
small quark masses, we nevertheless review here the theoretical considerations relevant 
for that situation, given that reaching this regime plays a central role for the general 
strategy presented in this paper. 

It was realised by Gasser and Leutwyler [42] (see also [43]) that, in a finite volume, 
low-momentum modes become increasingly important as the quark masses approach 
the chiral limit. Through a re-definition of the chiral counting rules one defines the 
so-called e-expansion, which represents a systematic low-energy description of QCD in 
a finite volume for arbitrarily small quark masses. Accordingly, the field U £ SU(4) is 
written as 

U(x) = e i2 ^ x)/F U , (5.11) 
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where £ describes non-zero momentum modes only, while Uq is a constant SU(4) matrix 
collecting the zero modes. The integration over Uq must be carried out exactly when 
mEF ;$ 0(1)> where m is the quark mass, while the integration over the non-zero 
modes may be carried out perturbatively, provided that F in units of the box size is 
large: FL ^> 1. The power counting rules for the e-expansion are then 

F~0(1), <9 M ~0(e), L„~0(l/e), £ ~ 0(e), m ~ 0(e 4 ), (5.12) 

where is the box length in direction /j>. Note that the quark mass counts as four 
powers of momenta, rather than two, as in standard ChPT in infinite volume. This 
has the important consequence that additional interaction terms, which appear beyond 
leading order, are suppressed if they contain one or more powers of the quark mass. 
For instance, if one works at next-to-leading order, including corrections of 0(e 2 ), the 
physical pion mass and decay constant, M n and F n , differ from their leading order 
values only by terms of relative order e 4 [44], owing to the fact that no higher-order 
terms contribute to the action at order e 2 . Furthermore, no additional interaction terms 
are generated in the effective weak Hamiltonian at order e 2 , and hence the knowledge 
of the associated LECs is not required in the analysis of the AI = 1/2-rule [45]. The 
fact that higher-order terms do not contribute at next-to-leading order in the e-regime 
represents an enormous simplification over the standard chiral expansion, where a large 
number of additional interaction terms arises at next-to-leading order in the effective 
weak Hamiltonian [24]. 

The LECs g^ 1 , g^ in the effective weak Hamiltonian can then be determined via 
the matching of correlation functions in the e-regime, i.e. in a finite volume, close to the 
chiral limit. Such a procedure has already been applied with some success for quantities 
such as the quark condensate S [46-48] and the pion decay constant [49,50]. 

5.3 Correlators in the e-regime 

We now proceed to define the relevant correlation functions. Details of the calculation 
are presented in Appendix[Dl We also refer the reader to ref. [45] , which gives a detailed 
account of the same calculation in the more conventional SU(3)l x SU(3)r case. 

We focus on correlation functions of two left-handed currents at Euclidean times 
xq and yo, respectively, and operators Of, Of at z = 0. Their definitions are 

Jd 3 x(j a (x)J \0)), (5.13) 

J d 3 x J d 3 y (j Q a (x)[dt(0)U^(y)) , (5.14) 

J d 3 x J d 3 y (j a (x)[df(0)U 7 sJ^y)) , (5.15) 

where the integrations are performed over the spatial volume L1L2L3. Since the pro- 
jected operators [0 1 ] a p-ys are linear combinations of [0i] a/ 3 7 <5, the correlation function 



C ab (x ) = 

[C 2 {xo,yo)] a ^5 = 
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a'p'yS ls obtained as a sum of correlators of the individual operators 0\- The same 
holds for [Cf]% lS . 

The general expressions for the correlation functions are presented in Appendix iDl 
and [45] . Here we only quote the results for a specific set of parameters, which corre- 
sponds to the unphysical situation of a light charm quark. We choose a diagonal mass 
matrix, M = diag(m, m, m, m), and work in a four-volume V = T x L 3 . Keeping in 
mind that = F and M n = (2mS/F 2 ) 1 / 2 at order e 2 , we find 

C a \x Q ) = Tr(T a T b )^ jl + [iV f ((3 ip - 3 / 2 - k 00 ) + uC e (u/2) h^Xo/T)] } , 

(5.16) 

where p = T/L, and Cg(u/2) is defined in eq. ()D.8j) . Here we keep explicit factors of 
Nf, as this will turn out to be useful when the quenched theory is considered later. The 
function hi and the "shape coefficients" /3i and feoo are specified in Appendix El The 
variable u is given by u = M^F^V = 2raYA/ [42,44,51]. For the same parameter choice, 
and assuming that the generators T a and T b are chosen according to eq. (|D.6f) . we find 
for Cf(x Q ,y Q ): 

^(xcyo) = ^{l + -^p 3 [2(iV f + (r)(Ap- 3/2 -feoo) (5-17) 

+ uC7 e (u/2)(/n(ar /r) + h^yo/T))] |, a = ± 

where we have suppressed flavour indices for brevity. In order to study the convergence 
properties at this order of the expansion, we plot Cf (xo,yo) i n Fig- D(left panel) for a 
volume with time extent T = 3.0 fm and two different spatial box sizes L, corresponding 
to 1.9 and 2.1 fm, respectively. We have also set 9 = here. One clearly sees that the 
corrections at order e 2 are quite large for L < 2.0 fm. 

The above results are obtained in a fixed ^-vacuum. By performing a Fourier 
transform in 9 one can obtain averages in sectors of "fixed topology", characterised by 
an index v. In ref. [52] the special role of topology in the e-regime was emphasised, 
noting that certain quantities depend quite strongly on the index near the chiral limit. 
Recent developments in the understanding of the role of topology show that correlation 
functions in fixed topological sectors can be defined also in QCD [53,54]. Although this 
result is universal [54], different topological sectors can be identified through the index 
of the fermion operator only if fermionic discretisations with an exact chiral symmetry 
are used. 

When considering averages in sectors of fixed v, poles in the quark mass m are 
expected to occur in certain correlation functions. Whether or not such poles also 
appear in quantities considered here is what we discuss in the following. Note that in 
the regime of larger quark masses (conventionally called p-regime), topology does not 
play a major role, and it is commonplace to present predictions of ChPT for fixed 9 
rather than v. 
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Figure 1: Left: the correlator Ci(xo,yo) of eq. 15.1711 . normalised to the tree-level value (solid line). 
The time extent is T — 3.0 fm, the quark mass m = 5MeV, the ^-coordinate Xq/T — —1/4, the pion 
decay constant F„ = 93MeV, and the chiral condensate E = (250 MeV) 3 . Right: the same observable 
with the spatial extent L — 2.1 fm, but in an ensemble with a fixed topological charge v. 



The transition to sectors of fixed topology simply amounts to substituting C$(u/2) 
by C v (u/2), whose definition is specified in AppendixO As can be seen from Fig.^right 
panel) the time behaviour of the correlator is indeed strongly modified in the small-mass 
region. Although the correlators do not develop any poles for m — > (because C u (u/2) 
is multiplied by u), their time dependence does not vanish in this limit, for v ^ 0. 

So far we have ignored the correlation function [C^ 1 ]. The reason is that it arises 
only at order e 8 , owing to the fact that it contains two powers of the mass matrix 
each contributing a factor of order e 4 . Actually, if a diagonal mass matrix and flavour 
assignments as in the derivation of eq. (|5.17|) are chosen, [C^ 1 ] vanishes identically. In 
order to isolate the "pure QCD" contribution to the AI = 1/2 enhancement at next-to- 
leading order, it is thus sufficient to concentrate on the LECs g x , dropping the operator 
[O2 ] a pj6 altogether. 

The determination of gf and is greatly facilitated by considering suitable ratios 
of correlation functions. If we define 

(xo ' yo) = Wo)^w (5 - 18) 

choose T flt = (T a )^, T fet = (T b )\ with T a , T b as in eq. (|D.6|) . and insert the expressions 
derived above we obtain 

K a {x , ya) = 1 + j^p 3 {/?i/T 3/2 - k 00 } , a = ±. (5.19) 
Thus, at order e 2 the dependence on m, xq, yo and — if the correlations are computed 
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for fixed topology — also on v drops out in the ratios K^. If we write 



K a (x ,y ) = l + aR(x , y ), * = ±, (5.20) 

we can plot the expression for R(xo,yo) (which is independent of xo,yo at order e 2 ) as 
a function of the volume, see Fig. [21 Another important ratio is 

H{XQ,y ) = -~7T. ~ = TTT7 V 5 ' 21 

Cf{x ,y ) K+{x ,y ) 
In terms of the one-loop correction R(xo,yo) it is given at 0(e 2 ) by 

H(x , y ) = l- 2R(x , y ) = l- ~^p 3 {(3i P ~ 3/2 - k 00 } , (5.22) 

and the right panel in Fig. [21 shows a plot of H as a function of the volume for different 
geometries. The computation of the corresponding ratio of three-point functions in 
lattice QCD in the e-regime is then directly proportional to g± /gf: 

9i _ k i Ci(x ,y ) 1 

gf kf Cf{xQ,yo)H(xQ,y Q y 

where 

C±(x ,y ) = ^(jo(x)Qf(0) J (y)) (5.24) 

is the three-point correlation function in QCD. Here it is assumed that the Wilson 
coefficients kf' are known for the particular scheme in which the operators Q 1 are 
renormalised. 

Let us stress that, in the main result of this section, eq. (|5.22j) . the Goldstone 
zero-momentum mode contribution, Cg(u/2) or C v [u/2), has cancelled completely at 
this order. This means that the result in the e-regime at NLO can be obtained from a 
NLO calculation of the same quantity in the p-regime, by taking just the leading term 
in a Taylor expansion in the pseudoscalar mass nip = 2mS/F 2 . Note that, as long 
as the volume is kept fixed, the contributions from non-zero modes are infrared-safe, 
since the chiral logarithms encountered in infinite volume are replaced by logarithms 
of the volume. The fact that the expression for the ratio /Cf in the p-regime can 
be Taylor-expanded in mp for fixed volume can also be used in practice to obtain the 
result in the e-regime indirectly from numerical simulations in the p-regime, by fitting 
the numerically determined ratio C{~ /Cf to a constant (which is the result in the e- 
regime) plus 0(mp) corrections with some unknown coefficient. 

5.4 The quenched approximation 

Currently, numerical simulations of lattice QCD with Ginsparg- Wilson fermions are 
mostly restricted to the quenched approximation, owing to the large computational 
cost. Here we discuss the implications of quenching for our strategy. 
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Figure 2: The quantities R(xo,yo) (left panel) and H(xo,yo) (right panel) plotted as functions of the 
box size and geometry. For F we have here chosen a value indicated by quenched simulations at modest 
volumes, F ^ 103 MeV [50], but the dependence of R on F is only through an overall factor 1/(FT) 2 
and the curves have thus the same form for any other value as well. 



Quenched QCD can be treated within the framework of an effective Lagrangian, 
with the resulting low-energy expansion referred to as quenched Chiral Perturbation 
Theory (qChPT) [56,57]. The theoretical status of qChPT is, however, questionable. 
This manifests itself in the occurrence of infrared divergences in certain correlation func- 
tions. These divergences reflect, at least partially, the sickness of quenched QCD. Here 
we adopt the pragmatic assumption that — despite the fact that it is not an asymp- 
totic expansion of quenched QCD (for a fixed number of colours N c ) — quenched ChPT 
does describe the low-energy regime of quenched QCD in certain ranges of kinematical 
scales, where correlation functions can be parameterised in terms of effective coupling 
constants, the latter being defined as the couplings which appear in the Lagrangian of 
the effective theory. Under this assumption we now discuss the determination of the 
counterparts of gf in the quenched theory, via the matching of correlation functions 
computed in quenched QCD and qChPT. 

The main difference between ChPT considered in the quenched and unquenched 
theories is that in the former there is no decoupling of flavour singlets from the pseudo- 
Goldstone bosons. The dynamics of the singlet field <3>o must then be incorporated into 
the effective chiral Lagrangian, which requires additional LECs associated with the new 
interaction terms. To be more precise, we consider the quenched effective Lagrangian 
in the so-called "replica formalism" [58, 59] 

- ±£Tr \y e UM ] + MU^Uq 1 



^qucn 
E 



±F 2 Tr 



(0„17)V 
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In the replica method one introduces N v valence quarks that are considered indepen- 
dently from the number of sea quarks, Nf, which will eventually be taken to zero. In the 
quenched case the fields U are promoted from elements of SU(iVf) to those of U(iVf). 
The factor Ug is given by Ug = eyip(i0I]y v /N v ), where In v is the unit matrix in the 
subspace corresponding to valence quarks. In addition to the LECs F and S (which, of 
course, may differ from those in the full theory), there are new parameters, a and mo, 
which are associated with the singlet field <I>o- 

We have noted above that poles in the quark mass may appear in fermion propaga- 
tors if the theory is considered for fixed non-trivial topology. These poles are expected 
to be the same in the full and quenched cases. For fixed topology, the counting rules of 
the e-expansion remain unchanged when passing to the quenched theory. 

When the AS = 1 weak interactions are incorporated into the quenched setting, 
the larger symmetry group and the presence of the singlet field $o may allow for ad- 
ditional interactions that need to be taken into account at a given order of the chiral 
expansion. Without going into further detail, we note that the current J® and the 
operators 0\ and O2 need not be modified in the quenched case, provided that an addi- 
tional expansion in 1/N C is assumed, which is needed in order to justify the truncation 
in operators involving higher powers of <I>o- This truncation has also been invoked to 
arrive at eq. 1)5.25(1 . Furthermore, the effective weak Hamiltonian does not have to be 
supplemented by additional operators. The classification of operators according to the 
valence symmetry does not give rise to ambiguities like those discussed in [60] for op- 
erators that appear in the conventional SU(3)-case. The argument which proves the 
absence of these ambiguities is analogous to the one discussed in [45] for the 27-plet 
operator in the SU(3) case. We refer the reader to this reference for more details. 

We now give the expressions for the quenched analogues of the correlation functions 
in eqs. (|5.16|) and (|5.17j) for fixed topology. They are formally obtained by taking Nf — > 
in the general expressions for the correlators obtained from the Lagrangian in eq. (|5.25|) . 
For the two-point function this yields [45,61] 

C a6;qucn (*o) = Tr(T a T b )^ {l + -JL_p 3 [uCr n (u/2) h^/T)] } , (5.26) 

where 

a qu ° n = ^(ReTr v C/o), iC/o (5.27) 

and Tr v denotes the trace over the valence subgroup. The precise meaning of the 
average (ReTr v f7o)v,j7 is specified in [45]. Note that there is no dependence on the 
singlet couplings m(j and a. As a result, it is easy to see that the quenched results can 
be simply obtained from the unquenched expressions in the limit Nf — ► 0. In a similar 
manner one finds the result for the three-point function 

Cr qucn (*o,2/o) = ^{l + ^P 3 [MAp- 3/2 -A;oo) (5.28) 
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+ 



uCr>/2)(M*o/r) + hiyo/T)] 



One can easily see that the ratios K ± (xo,yo) an d H(xo,yo) defined above remain un- 
changed, and thus the relation in eq. (|5.23|) between the ratio g^ /gf and the correlation 
functions carries over to the quenched case without modification, at this order. The same 
is true for the expansion of the ratio Cj~ /Cf in nip at fixed volume, which was discussed 
at the end of Sect. 15.31 

6 Decoupling of the charm quark in Chiral Perturbation Theory 

We will now depart from the unphysical situation of a charm quark degenerate with the 
other light quarks and investigate its effects in the framework of ChPT. In particular, 
we will show how the AI = 1/2 and 3/2 amplitudes evolve as m c is increased, but stays 
below the QCD scale, in order for ChPT to remain applicable. The aim is to find a 
relation between the LECs gf 1 and the corresponding couplings in the theory where the 
charm quark has been integrated out. Here we only sketch the procedure and state the 
main results: a detailed account of the calculation is given in a separate publication [62]. 

6.1 Basic setup 

The effective Lagrangian £e of QCD with an SU(4)l X SU(4)r chiral symmetry was 
defined to leading order ineq. ()5.1|) . Since the decoupling of the charm quark is insen- 
sitive to infrared physics, we will, for simplicity, consider the theory in infinite volume. 
Unlike the situation encountered in the e-regime, counterterms must then be added 
to the effective Lagrangian, if higher orders in the expansion are considered. For our 
purposes it is sufficient to include the counterterm 

5£ E = L 4 Tr[d^[/t]Tr[ X t[/ + tf X ], (6-1) 

where the vacuum angle 9 has been set equal to zero and x ls defined as 

2SM , . 

X = ~YT- (6-2) 

Phenomenological estimates for the value of L4 are available in the physical SU(3)l x 
SU(3)R-symmetric case [63,64]. 

Weak interactions are incorporated through the effective Hamiltonian, which, at 
leading order, contains the operators 0\ and O2 given in eqs. (|5.6|) and ()5.7|) . We now 
consider a mass matrix of the form 

M = diag(m u , m d , m s ,m c ) (6.3) 

with m c 3> m u = = m s , noting that the restriction to degenerate non-charm quarks 
is irrelevant for our purposes. Furthermore we define 

= 2m u T, _ 2m c S 

Xuu — p2 > Xcc — p2 ' I"- 4 ] 
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If momenta p with p 2 > Xcc are excluded, then one expects the physics of the SU(4)- 
symmetric theory to be described by another effective Lagrangian from which the heavy 
scale has been integrated out, and which possesses an SU(3)l x SU(3)r chiral symmetry. 
Our task is now to derive the effective weak Hamiltonian of such a theory, assuming 
that the LECs gf 1 are already known in the SU(4)-symmetric case. In order to define 
our power-counting rules we assume the following hierarchies: 

Xuu < Xcc < (4ttF) 2 . (6.5) 

In accordance with these relations, we work at order Xcc/i^F) 2 in the chiral expansion, 
while dropping all terms of order Xuu/Xcc and Xuu/(^F) 2 . 

The form of the strong interaction part of the effective SU(3)l x SU(3)r theory 
is identical to C~e, but formulated in terms of SU(3) matrices and modified LECs. In 
order to distinguish them from their counterparts in the SU(4)l x SU(4)pj,-symmetric 
theory, we will denote them with a bar. The effective action then reads 

C E = \F 2 Tr Ud^df.W] - ±E Tr \pM ] + MW] , (6.6) 
where U £ SU(3) and 

M = diag(m u ,m d ,m s ). (6.7) 

Flavour indices are denoted by a, f3, .. . and take their values from the set (u, d, s). The 
construction of non-singlet weak operators transforming under 3* (g> 3* (g> 3 <S> 3 of SU(3)l 
proceeds along the same lines as in the SU(4)-case. By considering the power-counting 
rules introduced above, one finds that the SU(3)-counterpart of O2 is of relative order 
(Xuu/Xcc) 2 and can thus be dropped. This leaves 

[<5xW- - \F^(ud^)_ & (ud^) w (6.8) 
[0 3 ]^ = \F 2 t{UM^ + MU% a , (6.9) 

and the effective weak Hamiltonian, denoted by 7Y^ hPT , is a combination of these two op- 
erators. The matching proceeds by enforcing equality between correlation functions in- 
volving left-handed flavour currents. In analogy with the procedure outlined in Sect. 15. ll 
the latter are derived by defining covariant derivatives according to 

d^U - D„U = (<9 M + iA*f a )U, dJJ -> D^U = {dp + iA%f a )U, (6.10) 

where T a are the Hermitian generators of SU(3). 5 The left-handed currents J® and 
in the SU(4)- and SU(3)-symmetric theories are then obtained by taking functional 
derivatives with respect to A^. 

4 In order to make contact with the conventions of ref. [45], we note that there the operator Oi is 
denoted by O m , while our [(5;i] sd corresponds to Og of [45]. 

Or, more precisely, in the case of SU(4), Hermitian generators in the sub-algebra that generates 
SU(3). 
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In order to perform the matching to leading order in the weak Hamiltonian we 
require 



JZ(*)JZ(v)) = (JSWJto))^, (6-n) 



SU(4) \ v ' " W/ /SU(3) 

/ChPT/^N <rb(„.\\ _ / -fa/„A-T>ChPT 



^(*)W^W>)) gu(4) = {Jfa)n^ l (z)JZ(y)) svi3) , (6.12) 

where the expectation values are evaluated using the strangeness conserving Lagrangian, 
and it is understood that space-time separations are large compared with the scales set 
by Xcc and 4irF. The objective is then to compute the observables on the left-hand 
sides of eqs. (|6.11(l and 1)6.12(1 to the order in Xcc defined previously and find out how 
the LECs and operators on the right-hand sides must be adjusted. 

6.2 Matching of the SU(3) and SU(4)-symmetric theories 

By enforcing the matching condition on current-current correlators one obtains the 
relations between the low-energy parameters F and F, as well as £ and £. Here we 
omit all the details of the calculation and refer to [62]. At one-loop level, the pion decay 
constant in the chiral limit in the SU(3)-symmetric theory, F, is related to F via 



F 2 = F 2 i 1 - - 



1 Xc 



l n Z££_l_l 6 ( 4vr )2 L MSjl /g 13) 
2u z J J 



2 (4vrF) 2 L 2/2 

Here we have adopted the conventions of the MS-scheme, with /2 denoting the subtrac- 
tion scale. A similar relation can be derived for the parameters £ and £ [62]. 

When matching the effective weak Hamiltonians according to eq. (|6.12jl . it is con- 
venient to replace W^ hPT (2:) by the operator 0\ in the SU(4)l X SU(4)r theory and 
work out the set of operators in the SU(3)l X SU(3)r theory it leads to. Skipping over 
the details of the calculation, we simply quote the result: when moving from 7^ hPT to 
W^ hPT , the operator [Ox] must be replaced by 



+d 2 (6 ac 6 Sc [0 3 ] lh + S Pc 8 lc [6 3 ] & s) , (6.14) 
where the coefficients c\ , C2 and d 2 are given by 



Cl = l + 2K^, (6.15) 

3 Xcc , j 1 Xcc , /„ , „n 

C2 = -- A (A ^ In , d2 = -- In , 6.16 

4(4vrF) J Xcc 2(4vrF) 2 Xcc 

and SU(3) singlet structures have been omitted in eq. (|6.14|) . Here A x is some physical 
scale that will in general be different for c 2 and d 2 , since it incorporates the effect of 
the finite corrections (similar to the terms proportional to in eq. (|6.13|) . but for 
the weak interaction part [24]) that have not been included in this case. On the other 
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hand, K denotes a particular such coupling that could contribute at this order through 
the higher-dimensional operator 

^[Ox^TrtfU + rfx\- (6-17) 

Since this operator is of order Xcc/F 2 it must be included according to our power 
counting rules. Thus we find that C2 and c?2 contain a logarithmic enhancement in m c , 
while ci does not. 

A weak operator [Cy^^j of the form shown ineq. (|5.7|) can also contribute to the 
effective weak Hamiltonian 7Y^ hPT . For our choice of mass matrix and the adopted 
power counting scheme, there is indeed a tree-level effect of the form 

2 

- ^(S ac 6 Sc [6 3 ]^ + 6 Pc S lc [d 3 y), (6.18) 

which is of order m^. In the region of small m c this is parametrically smaller than the 
term with the same structure in eq. (|6.14|) . which is multiplied by d%. As the latter is 
formally of order m c , we drop the above contribution from now on. 

We are now in a position to work out the effective weak Hamiltonian in the 
SU(3)l x SU(3)r theory. To this end we start with the expression for 7i^ hPT in the 
SU(4)-symmetric theory, eq. (|5.8|) . drop the contribution from O2 for the reason just 
mentioned, and write the remaining terms in the form 

g 2 

[OlUpr, (6.19) 

where the projectors Pf and P% are defined in eqs. (|A.3j) and (|A.4|) . From the above 
discussion we know that, to leading order in m c , the matching is achieved by replacing 
0\ by the right-hand side of eq. 1)6.14]) . with the coefficients c\, C2 and cfo specified as 
above. The expression for 7i^ hPT is obtained by inserting eq. (|6.14|) into eq. (|6.19|) and 
performing the decomposition into irreducible representations of SU(3). 6 The resulting 
expression for 7Y^ hPT is rather lengthy, and here we simply quote the result for the 
physical choice of flavours, characterised by c suu d = 1 and c SCC( i = — 1: 



-% (9t+9T)l&zU}- (6.20) 
In this expression the operators [O l ] suu d and [T^/ - ]^ are defined as 

[0\ jsuud = \ {\P^\suud + ~ ^[OljsKfid) = f ([Ol]sndu + §[0l]suud) ,(6.21) 

[UtU = \[Oi]sRRd. (6.22) 
6 Our conventions for the SU(3) classification are listed in Appendix A of [45]. 



21 



The operator [O x } suu d transforms under the irreducible representation of dimension 27 
of SU(3), and in eq. (|6,21|) we have listed two equivalent forms which are found in the 
literature. The operators [72. J" and [03] s <2 both transform under the 8-dimensional 
irreducible representation. 

After inserting the expressions for the coefficients ci, C2 and d<i while keeping only 
logarithmic terms we finally arrive at 
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This result illustrates the effects of the decoupling of the charm quark. One observes 
that the coefficients of the octet part (i.e. the last two lines) contain an extra logarithmic 
factor, Xccln A 2 /xco which for Xcc A 2 is parametrically larger than the coefficient 
of the 27-plet (first line), even if the linear term of eq. (|6.15|) were kept in the latter. 
It is well known that the octet part only mediates AI = 1/2 transitions, whereas the 
27-plet contains both AI = 3/2 and 1/2 contributions. Thus, the logarithmic terms in 
front of the octet operators describe the departure of the AI = 1/2 amplitude from the 
SU(4)-symmetric situation. 7 

We stress that these findings apply to the case of a moderately heavy charm quark, 
which must be light enough for ChPT to be valid. Whether or not the observed en- 
hancement survives if the charm quark is tuned towards its physical value must be 
studied in a lattice simulation. 



7 Numerical results in the SU (4) -symmetric theory 

We now describe the findings of our numerical investigations in the theory with a light, 
degenerate charm quark, i.e. m u = rrid = m s = m c . Thereby we will gain information 
as to what extent physics at the intrinsic QCD scale of a few hundred MeV contributes 
to the enhancement of the AI = 1/2 amplitude. The effects of a heavier charm quark 
will be the subject of future investigations. 

In our numerical work we have used the quenched approximation, whose deficits in 
this context have been discussed in Sect. [3 Although one may call into question any 
physical interpretation of quenched numerical data, we regard the work presented here 
primarily as a study to yield valuable information for future simulations. 

7 Note, however, that the last term in eq. 16.231 does not contribute to physical kaon decays [1,65]. 
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7.1 General definitions 



Our task is the determination of the low-energy constants gf and g~[ . Of particular 
interest is the ratio <7i~/Si > which, according to eq. (|5.1()[) . determines the ratio of 
amplitudes Aq/A2 at leading order in ChPT. The determination of LECs proceeds as 
usual, by computing suitable correlation functions in simulations of lattice QCD and 
fitting the results to the expressions obtained in ChPT. At non-vanishing lattice spacing 
such a procedure is justified only if a fermionic discretisation is chosen which preserves 
chiral symmetry. Here we have used the Neuberger-Dirac operator, which was already 
discussed briefly in Sect. 13.11 

From the previous sections of this paper it is clear that we are particularly inter- 
ested in correlation functions of left-handed currents and four-quark operators. More 
specifically, we define the non-singlet left-handed current as 

[J»ix)\aP = & a ^P-^)(x), (7.1) 

where a, (3 denote generic flavour indices, and the modified quark field V is given in 
eq. 1)3. 5|) . Our earlier considerations show that we can concentrate on the operator Qf, 
whose definition we repeat here: 8 

Qf = {(a7 M P_t2)(u 7At P_d) ± (s 7M P_d)(^P_{t)} - (u^c). (7.2) 

We then consider the following two-point and three-point correlators: 9 

C(x ) = E<[ J o(*)W[</o(0)W, (7.3) 

X 

Cf(x ,y ) = J2({ J o(x)]du[Qf(0)][My)]u S ). (7.4) 

After performing the Wick contractions, the correlators can be expressed in terms of 
the quark propagator S(x,y), 

S a (x,y) = (&(a# a (iO) = {(1 - \aD)D^ a ) (x,y), (7.5) 

where 

Dm a = (1 - \am a )D + m a (7.6) 

is the massive Neuberger-Dirac operator with bare mass m a . The left-handed propaga- 
tor is then given by [66] 

PS a (x,y)P + = {l-\am a r 1 {P-( ri sD ma )- X P+}{x,y)- (7.7) 

8 Note that we have suppressed the superscript "bare" in comparison with eq. 13.711 . 
9 No summation over a, f3 is implied. 
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(a) (b) 

Figure 3: (a): the "Figure-8" diagram; (b): the "Eye" diagram. A grey circle denotes the insertion of 
. The left-handed current is inserted at the black circles. 

Since here we restrict ourselves to studying the case where all quark masses are degen- 
erate, we will omit from now on the flavour label on propagators. The expression for 
the two-point function C(xq) then reads 

C(x ) = J2(Tr{ l0 P-S(x,0^P +l0 P-S(x,0)P + }). (7.8) 

X 

Turning now to the three-point correlator, we show in Fig.|3]the two types of diagrams 
that are obtained after performing the Wick contractions. An important consequence of 
restricting to the SU(4)-symmetric case is the fact that contributions from the so-called 
"Eye" -diagram vanish identically. Thus, only diagrams of the "Figure-8" type must be 
considered, and the expression for Cf 1 in terms of quark propagators becomes 

3 

Ct(x ,yo) = EE{(Tr{7 / ,P_5(x,0)t 7o P_5(x,0)}Tr{ 7/ ,P_5(y,0)t 7o P_5(y,0)}) 

x,y (J.=0 

T (Tr{ 7M P_S(x,0)t 7o P_5(x,0) 7M P_5(y,0)t 7o P_5(y,0)})}. (7.9) 

Thus, correlation functions of operators which transform under irreducible representa- 
tions of dimensions 84 and 20 are obtained by taking appropriate linear combinations 
of colour-connected and colour-disconnected contractions. As an aside we remark that 
the correlator (xo,yo) can also be used to compute the S-parameter Bk- 

In accordance with the relation between Qf and 0\ (see eq. (|A.9|1 ). we note that 
the combination k x C-y (xq, yo) corresponds to 2g^Cf{xQ, yo) in the low-energy descrip- 
tion, while kfC^(xQ,yo)/[C(xo)C(yo)] corresponds to gfK ± (xo, yo) (see eqs. (j5.14j) and 

(EH)- 

7.2 Technical details 

In order to evaluate the correlation functions defined above, we have computed left- 
handed quark propagators on quenched background gauge configurations, using the 
Neuberger-Dirac operator with s = 0.4 (see eq. 1)3.2(1 ) . We have employed the numeri- 
cal techniques described in ref. [66], including the approximation of the inverse square 
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Lattice 


P 


L/a 


T/a 




L[im] 


configs. 


A 


6.0 


12 


40 


7 


1.12 


751 


B 


5.8485 


12 


30 


5 


1.49 


638 



Table 1: Simulation parameters, including the spatial and temporal lattice sizes, the number of low 
modes of the Dirac operator which are treated exactly, and the number of gauge configurations. 

root in D by a minmax polynomial, the determination of the index v of a given gauge 
configuration via the counting of zero modes, and the computation of a number of low 
modes of P a D^ DP a , where P a projects onto the chirality sector without zero modes. 
The calculation of the left-handed propagator was accelerated using low-mode precon- 
ditioning as described in [66]. Furthermore, the calculation was arranged such that the 
necessary Conjugate Gradient inversions were always performed in the chirality sector 
without zero modes. 

In several recent publications [50, 67] it was reported that correlation functions 
of left-handed currents show strong statistical fluctuations if the quark mass becomes 
of order or smaller. The origin of these fluctuations could be traced to the 

low modes of the Dirac operator. For instance, if m — (SI/) -1 the contributions of 
a few eigenmodes to a given observables can be substantial, and the intrinsic space- 
time fluctuations of the eigenmodes then induce a large variance in the Monte Carlo 
estimate. In ref. [50] we proposed an exact technique which is able to reduce these 
fluctuations significantly, and which involves taking volume averages of the contributions 
of a certain number of low modes, n\ ow , to the respective correlators. This technique, 10 
dubbed "low-mode averaging" , allows for the computation of two-point functions with 
controlled statistical errors in all topological sectors if ni ow < 7, for quark masses around 

(EF)" 1 . We note, though, that our technique cannot cure the fluctuations caused 
by extremely small eigenvalues of D, which are expected to occur with a non-negligible 
probability if m <C (SI/) -1 . As already remarked in [50], the solution to this problem 
may require the incorporation of the low-mode contribution to a particular observable 
into the importance sampling process. 

For this study we extended the technique of low-mode averaging to the three-point 
functions corresponding to the diagrams in Fig.|3J It turned out, however, that low-mode 
averaging for three point functions requires a more careful tuning of its free parameters, 
in order to observe a significant reduction of statistical errors. Preliminary results of 
a detailed study have shown that masses m < (SF) -1 can be reached if the number 
of low modes treated exactly is increased by a factor 2 — 3 compared to the values in 
ref. [50]. 11 For the purpose of this paper we have, instead, focused our attention on 
moderately light quark masses, corresponding to the p-regime. The parameters of our 
simulations are listed in Tabled 

10 Similar methods were discussed in refs. [69] and [70]. 

1 We thank C. Pena and J. Wennekers for their contribution on this point. 
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Lattice am amp aFp C 1 /Cj~ 

A 0.030 0.2248(30) 0.03434(37) 2.46(40) 

0.040 0.2522(28) 0.03527(36) 1.94(25) 

0.050 0.2772(27) 0.03623(36) 1.67(18) 

0.060 0.3005(26) 0.03719(37) 1.52(15) 

0.070 0.3224(27) 0.03815(39) 1.41(12) 

B 0.040 0.2628(28) 0.04149(45) 2.03(25) 

0.053 0.2964(28) 0.04244(46) 1.75(17) 

0.066 0.3268(28) 0.04339(47) 1.57(12) 

0.078 0.3529(27) 0.04426(48) 1.45(10) 

0.092 0.3816(27) 0.04527(49) 1.35(8) 



Table 2: Results for pseudoscalar meson masses and the ratio of three-point functions C 1 /Cf at 
several values of the bare quark mass. 



7.3 Results 

In order to determine pseudoscalar masses and decay constants we fitted the two-point 
function C(xq) to the expression [50] 

ni \ i ^ cosh[(T/2 - |x |)m P ] 

C(s ) = - 2 m P F P 2ginh[rmp/2] • (7-10) 

For lattice A these fits were performed for timeslices in the range 7 < xo/a < 19, while 
for lattice B we used 7 < x^/a < 14. The results are listed in Table |2J We note that 
the values obtained on lattice A for bare masses am = 0.04 and 0.06 can be directly 
compared with ref. [50], where a larger spatial volume corresponding to L = 1.49 fm 
was used. The results for the pseudoscalar decay constant agree within errors, while the 
pseudoscalar meson mass at am = 0.040 is larger by 4% (1.8a) on the smaller spatial 
volume, which indicates a small finite-volume effect in the pseudoscalar correlator. 

In order to isolate the asymptotic behaviour of the ratio of three-point functions, 
Ci(xo,yo)/C±(xo,yo), we first fixed xq at xq = Xfi x = T/4 (actually, Xfi x = 8a for 
lattice B) and looked for a plateau in yo around T — x^ x . Figure shows the quality 
of the plateau for both runs at two values of the bare quark mass. By fitting the ratio 
/Cq t° a constant for 28 < yo/a < 30 (lattice A) and 18 < yo/a < 25 (lattice B), 
respectively, we obtain the results listed in the last column of Table [2j 

Our results show that there is a clear numerical signal for the ratio of correlators 
proportional to <q~ /gj in the SU(4)l x SU(4)r symmetric theory. The typical statistical 
error in the studied range of quark masses is at the level of 10%, which should be 
sufficient for the purpose of establishing whether physics at the intrinsic QCD scale 
makes a significant contribution to the AI = 1/2 rule. 

In Fig. [3 we plot the results for the asymptotic value of C-f /Cf for both our lattices 
as a function of the bare quark mass in units of the hadronic radius tq [74]. At the 
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Figure 4: The ratio C 1 (a?fi x , yo)/ (xfl x , yo) as a function of yo /a for lattices A (left panel) and B (right 
panel). Open circles and full squares correspond to the smallest and biggest quark mass, respectively. 



smallest quark mass the ratio of correlators of operators transforming under irreducible 
representations of dimensions 20 and 84 is clearly greater than one. Furthermore, one 
observes a clear trend for this ratio to increase as the chiral limit is approached. 

At leading order in ChPT the connection between C± /C^ and the ratio g± /gt of 
LECs is given by 

gT _ K C{(x ,yo) 

at kf Cf(x ,y ) ren ' 

where the subscript "ren" reminds us that the four-quark operators in the lattice- 
regularised theory must be related to a particular continuum scheme. Thus, in order 
to estimate /gf it suffices to multiply C{" /C± by the relevant short-distance factors. 
When interpreted in this way, our results indicate a significant enhancement of the 
AI = 1/2 amplitude from long-distance QCD effects alone. However, the prediction 
of leading-order ChPT in eq. (|7.11|) implies that there should be no significant mass 
dependence in the ratio, a behaviour which is not reflected in our results. 

Clearly, NLO effects must be taken into account. In the p-regime they depend 
on many unknown effective couplings, while, as we have shown in this paper, they 
can be computed in the e-regime without adding any extra parameters (see eqs. (|5.23f) 
and 1)5.22(1 '). Since we lack any data in the e-regime, the NLO matching of C-f jC\ to 
g\ I g\ cannot be performed here. However, in order to illustrate how the NLO formulae 
might be used, we note that for lattice B the correction factor H(xo,yo) evaluates 
to 2.262, independent of xq and yo- Thus, if long-distance QCD effects alone were to 
enhance g^ /gf significantly, say, by a factor 2, then the ratio C-f jC^ would have to be 
as large as 4 — 5 in the e-regime. 

We note, though, that the geometrical factor p = T/L enters H(xo,yo), such that 
for space-time geometries with T ^> L the NLO corrections can become large. For 
instance, on lattice B we have p = 2.5, and one might expect that higher orders in the 
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Figure 5: The fitted value of the ratio C-j" /C^~ plotted versus the bare quark mass in units of ro- The 
symbols are slightly displaced for clarity. 

e-expansion cannot be neglected in this case (cf. Figs.n an d|2J). 

Therefore, while we observe a good numerical signal for the ratio Cf/Cj 4- , a quan- 
titative interpretation of our results is not yet possible, since our masses are not in a 
regime where computable NLO corrections can be applied. Future simulations should 
therefore concentrate on the determination of Cf jCy for smaller quark masses. More- 
over, choosing a small value of T/L, while keeping the spatial box length reasonably 
large (e.g. L ^ 1.5 fm) will lead to better convergence properties of ChPT in the e- 
regime. 

8 Summary and conclusions 

We have outlined a computational strategy which allows to disentangle and quantify 
several possible origins of the AI = 1/2 rule using lattice simulations of QCD. The 
main idea is to consider the AS = 1 weak interactions with an active charm quark and 
to compute correlation functions for K — > tt transitions, whose values are subsequently 
matched to the LECs which parameterise the weak Hamiltonian in the effective low- 
energy theory. While the mass-degenerate case, m u = = m s = m c , allows to 
quantify contributions to the AI = 1/2 enhancement from long-distance effects that are 
not due to a large charm quark mass, the effects of physics at scales around m c can be 
estimated by increasing m c above m u and evaluating the "Eye" -diagram in addition to 
the "Figure-8" -graph. From the behaviour of the ratio of amplitudes as a function of 
m c one will be able to infer whether QCD long-distance effects are merely amplified by 
the splitting between up and charm quark masses, or whether the AI = 1/2 rule is due 
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to modes at a specific energy scale, which should manifest itself as a sharp change in 
magnitude. 

A key ingredient in our strategy is the use of fermionic discretisations which preserve 
chiral symmetry at non-zero lattice spacing (Ginsparg- Wilson fermions). In this case 
the renormalisation patterns of the relevant dimension-6 operators which are found 
in the continuum carry over to the lattice theory without modification. Moreover, 
the matching to the expressions of ChPT can be performed in a conceptually clean 
manner. Furthermore, the e-regime of ChPT offers a firm theoretical framework for the 
determination of LECs: owing to the specific chiral counting rules of the e-expansion, 
no additional interaction terms are generated in the effective weak part at NLO. By 
contrast, for the conventional chiral expansion (p-expansion) the matching to QCD at 
NLO involves coupling terms whose coefficients are unknown. 

Our studies have revealed, though, that fairly large volumes are required in order 
to guarantee reliable results. The slow convergence of ChPT in the e-regime is indi- 
cated in Figs, n an d El which demonstrate that corrections are relatively small only 
for L > 2.0 fm. The use of asymmetric lattices slows down the convergence further. 
Since lattice simulations are typically performed with T > L, in order to be able to 
isolate the asymptotic behaviour of Euclidean correlation functions, whilst keeping the 
total number of lattice sites at a manageable level, this may render future numerical 
simulations relatively expensive. 

We have performed numerical simulations in the SU(4)-symmetric case in the 
quenched approximation, with masses corresponding to the p-regime. The results show 
that the matching of C{~ /C± to gi / gf at leading order is clearly insufficient. On the 
other hand, for the reasons explained above, adjusting the kinematical variables so that 
NLO matching can be performed reliably in the e-regime is numerically expensive. In 
fact, the computational cost of simulating the kinematical regime where the p-expansion 
at leading order is valid might be even of comparable size. 12 Which of these alternatives 
is to be preferred must be decided by future simulations. 

We have also investigated the effects of a heavier charm quark in ChPT. These 
studies have shown that there is indeed a larger contribution to AI = 1/2 transitions 
when the charm quark mass is increased above the masses of the u, d and s quarks but 
kept below the chiral scale A x . Since only moderately heavy charm quark masses can 
be studied safely in ChPT, numerical simulations must be performed to confirm this 
result. This will be the subject of our future simulations, in which we will evaluate the 
necessary "Eye" -diagrams. Furthermore, we intend to improve on the systematics of 
our earlier runs, by using lattices with smaller T/L and going closer to the chiral limit. 

12 Various strategies for matching simultaneously for the LO and for a number of NLO couplings in 
the conventional SU(3)-case have been outlined in [75]. 
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A Four-quark representations 



In this appendix, we provide for completeness some aspects of the SU(4) classification 
of four-quark operators. For further details see, for instance, the complete presentation 
in [29]. We consider an operator O a/ g 7 5 which transforms under the 4* <g> 4* ® 4 ® 4 
representation of SU(4) and decompose it into its irreducible parts. To this end we 
define the projected operators 

Oaf3y6 = V~'\)oiP'y5;\vptO\vpt-i (A-l) 

where a = ±1. The projectors P[ , P% are given by 

{Pl)a(3~iS;\vpT = i(Sa\Sp u + Cr5 av 5(3x){b-ypb&T + Cr5j T 5sp), (A. 3) 

a/3"/S;Xupr = $a\$ f3u$~fp$8T + (4 _|_ 2cx)(4 + d) ^^^^ ^^aS^fj^^Xp^ur 

(6 a7 8i3 v 5 ST dxp + S l 3s5 a x5 lp 5 UT + oS aS Sp v S 7T Sxp + ^S^SaxSspS^r). (A.4) 



4 + 2o- 
Furthermore we define 

S* = OZ Xk x, (A.5) 

Ray = OaXyX ~ J^ayS" . (A. 6) 

Any four-quark operator a p 7 s can then be decomposed into irreducible representations 
according to 

+ 4 + 2er (Sa-yRps + SpsRay + cr6 a5 R^ + aS^R^g)^ . (A.7) 

The operators ^ i?„ 7 and transform according to representations of the follow- 
ing dimensions: 

°UrS ■ 84, R% : 15, 

<W ■ 20, S± : 1. (A - 8) 

In eq. I|2.7jl . the operators have been defined for the relevant physical flavour as- 
signments. Their relations to the projected operators Of 1 are given by 

Qf = 2([Of] suud - [Of} sccd ), (A.9) 

where Of are obtained from the generic four-quark operator [Oi] aj g 7< 5 via eqs. (|A.lj) 
and (|A~2|) . 
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B Exact chiral symmetry and CP 



In this appendix we report the argument spelled out in [31], which shows that in lattice 
QCD with Ginsparg- Wilson fermions, the CP symmetry remains a useful tool in the 
analysis of operator mixing, if a chiral operator basis is chosen. Assuming that D 
satisfies the Ginsparg- Wilson relation, we note that the quark action 

S F = a*J2fa x )l D ^( x ) +i>{x){MP + +M^P_)[{l-\aD)i>](x)) (B.l) 

X 

is invariant under the exact SU(4)l x SU(4)r symmetry group, provided the quark mass 
matrix M is transformed according to the 4 ® 4* representation. If the mass matrix 
M is chosen real and diagonal, the action is also invariant under parity and charge 
conjugation. The restriction to this case is not important, since chiral transformations 
can always be used to diagonalise M. However, four-quark operators such as 

0%h6 = Wa^P-^M^P-^) ±(l~6) (B.2) 

do not transform in a simple way under CP. At the level of correlation functions of local 
fields, this is not a problem, because one can substitute 

${x) = (1 - \-aMY 1 j^O*) " } (B.3) 

and perform a partial integration in the functional integral to eliminate the second term 
on the right-hand side. One then obtains the same correlation functions up to contact 
terms, but with the modified field ip replaced by (1 — ^oM)" 1 ^. 

In the case of the operators s one can show that all contact terms vanish after 
contracting the Dirac indices. When inserted in correlation functions of local fields at 
non-zero distances, one obtains the operator identity 

Otfrs = (! - ham^il - \arns)- 1 

x{(V^P_V 7 )(?/37^-V>5) ± (7 <- *)}■ (B.4) 
Applying a CP transformation yields 

± (X) CF : (1 - ^ m ^ 1 ~ jg^g) Q± (~) ( B5 ) 



where x = (xo, —xi, —X2, —xs), and it is understood that the rule only applies in local 
correlation functions. The argumentation in the case of two-quark operators proceeds 
along exactly the same lines. 
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C Perturbative results for the anomalous dimensions 7 <T (g) 

For convenience we recall that the RG /3-function has a perturbative expansion 

oo 

P(g) = -g 3 Y, b x9 2k , (ci) 

fc=0 

with the universal one- and two-loop coefficients: 

1 /ll 2 \ 

b = TT^ ^ N c ~ o ) , (C2) 



(4vr) 2 V 3 L 3 J 



1 /34, f2 „ f 13 



61 = (i^U^-My^-^jJ- (c - 3) 

The anomalous dimension of Qf has a perturbative expansion of the form 

oo 

f(9) = 5 2 E* a (C4) 

The one-loop coefficients 7q were first obtained by Gaillard and Lee [1], 

1 6a(N c -a) 

N c • (C - 5) 

The two-loop coefficients depend on the scheme and on the way 75 is handled. The first 
such two-loop computation was performed by Altarelli, Curci, Martinelli and Petrarca 
[33] using dimensional reduction [71]. The computation in the HV(MS) scheme was 
carried out much later with the result [34]: 



Tf _ 1 (N c -a) 



88 at2 01 16 AT AT f 157 AT 57 28 Ar 
■— iV c 2 + 21 + — N c Nf - a \ —N c + — - —N f 



(4vr) 4 2N C [ 3" c ' ~ ' 3 " I 3 " u ' N c ' 

Denoting the operators in the HV scheme by (Qf )hv, the renormalised operators of any 
other satisfactory scheme S should be related to these merely by a finite multiplicative 
renormalisation, 

(G?)s = *s ff G/)(Gi)Hv. (C7) 
Assuming that the normalisations are such that X§ = 1 at tree level, the coefficients 
have a perturbative expansion of the form 

00 

X^g) = l + J2^9 2k - (C8) 
fc=i 

The one-loop coefficients have been computed in various schemes, e.g. in the RI scheme 13 
one has 

x W = - { 4 ir) 2 { Nc 2 N c a) [4JVc-6q{l-41n(2)}] . (C.9) 



13 This is defined by demanding that the renormalised 4-point vertex function in the Landau gauge 
at equal external momenta p with p 2 = /i 2 (where /i is usually set equal to the standard renormalisation 
scale in the MS scheme) be equal to the tree level function. 
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For Ginsparg- Wilson fermions, as denned in Sect. one has 



*GW;1 = ( ^ )2 {NC 2Nc a) [XOV + X lNc ] , (CIO) 



with 



X = 4B y - 2B S + 2B^, + 6 , (C.ll) 
Xi = 2£ v + 2^-4, (C.12) 

where the functions By, 5g, B^ depend on the parameter s, and some numerical values 
can be found in Table 1 of ref. [39] . Corresponding values for the Xi are given in TableOl 
The anomalous dimensions of the operators in scheme S are related to those in the 
HV scheme by 

7s CT = 7 CT + /?(2)|^ln* s CT ( 5 ). (C.13) 

Denoting the coefficients in the perturbative expansion of 7g in the MS coupling (as in 
Eq. (|C.4|) ) by 7g. k , the lowest order coefficients are given by 

7s ; o = To, (C14) 



7S ;i = 7?-2& xg ;1 . (C.15) 

For completeness we mention that for the lattice computations it is often convenient 
to write the expressions in Sect. E] in terms of lattice couplings. In this case, as an 
intermediate step, it is useful to know the relation of the MS coupling to the bare 
lattice coupling: 

4^ a ~ 1 )=2o + ^/o 4 + ---- (C16) 

To obtain the anomalous dimensions expressed in terms of the lattice coupling to two- 
loop level, it is sufficient to know the coefficient d\ which takes the form: 

d x (t) = -8tt6 ln(i) + di + d u N f . (C.17) 

The term d\o depends on the pure gauge action; e.g. for the Wilson gauge action it was 
computed long ago by A. and P. Hasenfratz [72] and is given by 

dio = -^r + ^c, (C.18) 
fci = 2.135730074078457(2). (C.19) 

The term d%\ in l)C.17jl for GW fermions depends on the parameter s and is given by 

5 



du(s) = Air 



72vr 2 



(C.20) 



where numerical values for the functions dn(p),i = 1,2, with p = 1 + s, are given in 
Table 1 of ref. [73]. Corresponding values for du(s) are given in Table £3 
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s 


X 


Xi 


du 


-0.8 


-462.169 


-472.573 


0.287181 


-0.7 


-292.939 


-302.192 


0.256186 


-0.6 


-209.059 


-217.347 


0.235334 


-0.5 


-159.227 


-166.669 


0.219891 


-0.4 


-126.366 


-133.048 


0.207873 


-0.3 


-103.174 


-109.157 


0.198305 


-0.2 


-86.004 


-91.337 


0.190679 


-0.1 


-72.836 


-77.557 


0.184743 


0.0 


-62.460 


-66.599 


0.180413 


0.1 


-54.108 


-57.689 


0.177735 


0.2 


-47.270 


-50.313 


0.176883 


0.3 


-41.593 


-44.113 


0.178175 


0.4 


-36.827 


-38.836 


0.182126 


0.5 


-32.787 


-34.295 


0.189537 


0.6 


-29.338 


-30.351 


0.201685 


0.7 


-26.376 


-26.896 


0.220694 


0.8 


-23.820 


-23.848 


0.250421 



Table 3: Values of functions Xj appearing in eq. (|C.10|) . and d\\{s) appearing in 
eq. 



D Correlators in Chiral Perturbation Theory 

In this appendix we describe the calculation of correlation functions of left-handed 
charges and operators Of for a finite periodic box with volume V = T ■ L\ ■ L2 ■ L3, at 
order e 2 in ChPT. The same calculation in the SU(3)-symmetric case has already been 
published in ref. [45], which can be consulted for further details. In particular, we note 
that the sets of diagrams which appear at this order in the e-expansion, are identical to 
the SU(4)-case considered here, and are depicted in Figs. 1 and 3 of [45]. 
The result for the correlation function C ab at order e 2 is 



C 



ab 



u f 2 r N f 

Tr(T^)-[l + ^ 



Pi 



T 2 k, 



00 



yV2 



V 



+ 



(D.l) 



where Nf = 4, and the function hi is given by [55] 



hi{r) = \ (| 



I) 2 
2) 



j_ 

12 



(D.2) 



The definitions of the coefficients (5\ and /cooi which depend on the geometry of the box, 
can be found in [51]. Examples of numerical values are listed in Table H 
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T/L 


fk 


&O0 


32/32 


0.14046 


0.07023 = 


32/28 


0.13872 


0.07826 


32/24 


0.13215 


0.08186 


32/20 


0.11689 


0.08307 


32/16 


0.08360 


0.08331 


32/12 


0.00582 


0.08333 


32/8 


-0.21510 


0.08333 



Table 4: Some numerical values for /3i and fcoo, for geometries with L\ — L2 — L3 = L. 



For [Cj (xq, yo)]a« 7< 5 the general expression reads 



[Cf(x ,yo)}l% 5 



A 



+ 



ab,cr 



l + (AT f + (T )_ 



a/375 4 <p2 

^(ReTr[Mt/ e^I 
iV f F 2 \ 1 J /e,u 



T 2 k 



00 



V 



h 



(?) 



(D.3) 



Contact terms have been dropped in this expression, and we may assume that xq, yo 
are on opposite sides relative to the origin, i.e. xq < 0, yo > 0. The flavour structure is 
encoded in the tensor 



A 



ab,a 

a/375 



1 

2(4 + 2cr) 



7/3 5a 



(4 + a)(4 + 2a) 



+ a5 ia S Sf s)Tt(T a T b ) 



(<5 7/3 {T a , T b } 5Q! + Ssa{T a , T b } l/3 
+ a5 ia {T a ,T b } S(3 + a5 S p{T a ,T b } 



-ya j j 



(D.4) 



where 



T'farpb} rpa rpb I rpa rpb 



(D.5) 



1 7a- 1 5/3 7a 5/3 "i - 5/3 7a - 

In order to study the case of a light charm quark, one can choose a diagonal, degenerate 
mass matrix and generators such that 



rpa 
1 TUJ 



(D.6) 



where a, /3, 7, 8 are assumed all fixed and different. In this case one finds that A ab = 1/2, 
and for a symmetric spatial volume and vanishing vacuum angle 6 one recovers the 
expressions in eqs. ()5.16j) and (|5.17j) . 

Correlation functions in sectors of fixed topology, characterised by an index u, can 
be obtained from those for fixed vacuum angle via Fourier transform. If we assume a 
diagonal mass matrix and define 



Z e (u/2) 



3 (it/2)Re Tr[t/ exp(i0/N f )] 



1 



Uq 



2~ 



Z u (u/2) = --l ae e - ive Z e {u/2), (D.7) 

27T Jo 
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the quantity Cg(u/2) in eqs. (|5.16|) and (|5.17|) gets simply replaced by 

Ce{u/2) = ^-^lnZ (u/2) — > A|-ln^(u/2) = C v {u/2). (D.8) 
In the small and large mass regions, C u (u/2) is approximately given by 

CAu/2) « { 2 I"|/"' »«J- (D.9) 
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